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ABSTRACT 

We collect data at well sampled frequencies from the radio to the 7-ray range for the 
following three complete-samples of blazars: the Slew Survey and the 1 Jy samples of 
BL Lacs and the 2 Jy sample of Flat Spectrum Radio-Loud Quasars (FSRQs). The 
fraction of objects detected in 7-rays (E > 100 MeV) is ~ 17 %, 26 % and 40 % in the 
three samples respectively. Except for the Slew Survey sample, 7-ray detected sources 
do not differ either from other sources in each sample, nor from all the 7-ray detected 
sources, in terms of the distributions of redshift, radio and X-ray luminosities and of 
the broad band spectral indices (radio to optical and radio to X-ray). 

We compute average spectral energy distributions (SEDs) from radio to 7-rays for 
each complete sample and for groups of blazars binned according to radio luminosity, 
irrespective of the original classification as BL Lac or FSRQ. 

The resulting SEDs show a remarkable continuity in that: i) the first peak occurs in 
different frequency ranges for different samples/ luminosity classes, with most luminous 
sources peaking at lower frequencies; ii) the peak frequency of the 7-ray component 
correlates with the peak frequency of the lower energy one; iii) the luminosity ratio 
between the high and low frequency components increases with bolometric luminosity. 

The continuity of properties among different classes of sources and the system- 
atic trends of the SEDs as a function of luminosity favor a unified view of the blazar 
phenomenon: a single parameter, related to luminosity, seems to govern the physical 
properties and radiation mechanisms in the relativistic jets present in BL Lac objects 
as well as in FSRQ. The general implications of this unified scheme are discussed while 
a detailed theoretical analysis, based on fitting continuum models to the individual 
spectra of most 7-ray blazars, is presented in a separate paper (Ghisellini et al. 1998). 

Key words: quasars: general - BL Lacertae objects: general - X-rays: galaxies - 
X-rays: general - radiative mechanisms: non-thermal - surveys 



1 INTRODUCTION 

The discovery of BL Lac objects and the paradoxes associ- 
ated with their violent variability led to a major step for- 
ward in the theory of Active Galactic Nuclei (AGN) , that is 
to the concept of relativistic jets. Flat spectrum, radio-loud 
quasars (Angel & Stockman 1980) share basically all of the 
properties of BL Lac objects related to the presence of a 
strong non-thermal broad band continuum, except for the 
absence of broad emission lines. Hence the common desig- 
nation of blazars proposed by Ed Spiegel in 1978. 

It was initially supposed that BL Lacs represented the 
most extreme version of FSRQs, i.e. those with the most 



highly boosted continuum. Instead, it has been recognized 
later (e.g. Ghisellini, Madau & Persic 1987; Padovani 1992a; 
Ghisellini et al. 1993) that the amount of relativistic beam- 
ing and the intrinsic power in the lines are lower in BL 
Lacs than in FSRQs, implying some intrinsic difference be- 
tween the two classes. Differences are also found in the 
extended radio emission and jet structure (e.g. Padovani 
1992a; Gabuzda et al. 1992). Nevertheless the continuity 
of several observational properties including the luminos- 
ity functions (Maraschi & Rovetti 1994) , the radio to X-ray 
SEDs (Sambruna et al. 1996) and the luminosity of the lines 
(Scarpa & Falomo 1997) suggests that blazars can still be 
considered as a single family, where the physical processes 
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are essentially similar allowing for some scaling factor(s). 
The identification of these scaling factors would represent 
a substantial progress in the understanding of the blazar 
phenomenon. 

A special class of BL Lacs was found from identification 
of X-ray sources. X-ray selected BL Lacs (XBL) differ from 
the classical radio-selected BL Lacs (RBL) in a lesser degree 
of "activity" (including polarization), in the radio to optical 
emission and in the relative intensity of their X-ray and 
radio emission. This led to the suggestion that the X-ray 
radiation was less beamed than the radio one and that XBLs 
were observed at larger inclination to the jet axis (e.g. Urry 
& Padovani f995 for a review). 

Giommi & Padovani (1994) quantified the differences in 
SEDs between XBLs and RBLs, and Padovani & Giommi 
(1995) introduced the distinction between 'High-energy 
peak BL Lacs' (HBL) and 'Low-energy peak BL Lacs' 
(LBL), for objects which emit most of their synchrotron 
power at high (UV-soft-X) or low (far-IR, near-IR) fre- 
quencies respectively. Quantitatively a distinction can be 
done on the basis of the ratio between radio and X-ray fluxes 
(see also §3.2.2). We will use the broad band spectral index 
qrxQ and call HBL and LBL objects having qrxSO.75, 
>0.75, respectively. Giommi and Padovani also proposed 
that HBL represent a small fraction of the BL Lac popula- 
tion and are numerous in X-ray surveys only due to selec- 
tion effects. An alternative hypothesis (Fossati et al. 1997) 
relates the spectral properties to the source luminosity in 
such a way that low luminosity objects (with high space 
density) are HBLs while high luminosity objects (with low 
space density) are LBLs. 

We will include here X-ray selected BL Lacs together 
with "classical" BL Lacs in the blazar family, again assuming 
that the basic physical processes by which the continuum is 
produced are common to the whole family. 

The detection by EGRET, on board the Compton 
Gamma-Ray Observatory (CGRO), of many blazars at 7- 
ray energies (E > 30 MeV) revealed that a substantial frac- 
tion and in some cases the bulk of their power is emitted in 
this very high energy band. The 7-ray emission is therefore 
of fundamental importance in the SED of blazars. 

From the theoretical point of view the radio to UV 
continuum is universally attributed to synchrotron emis- 
sion from a relativistic jet, while a flat inverse Compton 
component due to upscattering of the low energy photons 
is expected to emerge at high energies as originally dis- 
cussed in Jones, O'Dell & Stein (1974). The latter pro- 
cess is therefore a plausible candidate to explain the 7-ray 
emission. The soft photons to be upscattered could be ei- 
ther the synchrotron photons themselves (synchrotron self- 
Compton process, SSC, e.g. Maraschi, Ghisellini & Celotti 
1992; Bloom & Marscher 1993) or photons produced by the 
disk and/or scattered /reprocessed in the broad line region 
(Blandford 1993; Dermer & Schlickeiser 1993; Sikora, Begel- 
man & Rees 1994; Ghisellini & Madau 1996). Understanding 
whether/how the 7-ray properties differ among subclasses 
is essential to assess the role of different mechanisms and to 



* Hereinafter we define spectral indices as F„ <x u~ 01 . In broad 
band indices radio, optical and X-ray fluxes are taken at 5 GHz, 
5500 A and 1 keV, respectively. 



verify whether the idea of blazars as a unitary class can be 
maintained. 

Here we study the systematics of the SEDs of blazars 
using data from the radio to the 7-ray band. We confirm 
and extend previous results of Maraschi et al. (1995) and 
Sambruna et al. (1996) by: i) extending the SED to the 7- 
ray range; ii) using a much larger complete sample of FSRQ; 
iii) using the richer and brighter sample of X-ray selected 
BL Lacs recently derived from the Slew Survey. We use the 
available 7-ray data for each sample but also indirect in- 
formation derived from the 7-ray detected (not complete) 
sample discussed by Comastri et al. (1997). Since we find 
that the continuity hypothesis among blazars holds we also 
consider a merged "global" sample subdivided in luminosity 
bins irrespective of the original classification of the objects. 

The structure of the paper is as follows. In section 2 
we describe how the data for the SEDs were collected and 
treated for each sample. The 7-ray properties of different 
samples are also discussed. In section 3 we build average 
SEDs for the three sub-samples and for the global sample 
subdivided according to luminosity and present our results. 
These are discussed in section 4, while our conclusions are 
drawn in section 5. 



2 THE DATA 
2.1 The samples 

We decided to consider the following three samples of 
blazars: the Slew Survey Sample and the 1 Jy sample of 
BL Lac objects and the FSRQ sample derived from the 2 Jy 
sample of Wall & Peacock (1985), motivated by the need of 
completeness, sufficient number of objects and observational 
coverage at other frequencies, as detailed below. 

2.1.1 BL Lacs, X-ray selected: the Slew survey sample 

The Einstein Slew survey (Elvis et al. 1992) was derived from 
data taken with the IPC while the telescope scanned the 
sky in between different pointings. It has limited sensitivity 
(flux limit of ~ 5 x 10~ 12 erg cm -2 s -1 in the IPC band 
[0.3 — 3.5 keV]), but covers a large fraction of the sky (~ 
36000 deg 2 ). In a restricted region of the sky Perlman etal. 
(1996a) selected a sample of 48 BL Lacs (40 HBL, 8 LBL) 
which can be regarded as being practically complete. This 
is the largest available X-ray selected sample of BL Lacs. 
The redshift is known for 41 out of the 48 objects, and 8/48 
have been detected at 7-ray energies, (6 with EGRET, E> 
100 MeV, 1 with Whipple, E> 0.3 TeV, 1652+398, and 1 
with both instruments, 1101+384). 

2.1.2 BL Lacs, radio selected: the 1 Jy sample 

This is the largest complete radio sample of BL Lacs com- 
piled so far. The complete 1 Jy BL Lac sample was derived 
from the catalog of extragalactic sources with Fsghz > 1 
Jy (Kiihr et al. 1981) with additional requirements on radio 
flatness (o?r < 0.5), optical brightness (my < 20) and the 
weakness of optical emission lines (EWa < 5 A, evaluated 
in the source rest frame) (Stickel etal. 1991). This yielded 
34 (2 HBL, 32 LBL) sources matching the criteria, 26 with a 
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redshift determination and 4 with a lower limit on it (Stickel, 
Meisenheimer & Kiihr 1994). Out of these 34 objects, 9 have 
been detected at 7-ray energies (8 with EGRET, plus 1 with 
Whipple, 1652+398). 

2.1.3 Flat Spectrum Quasars: Wall & Peacock sample 

For FSRQs we considered the sample drawn by Padovani & 
Urry (1992) from the "2 Jy sample" (Wall & Peacock 1985), 
a complete flux-limited catalogue selected at 2.7 GHz, cov- 
ering 9.81 sr, and including 233 sources with -F2.7GHZ > 2 
Jy, and Qr < 0.5. It consists of 50 sources with almost com- 
plete polarization data, of which 20 are detected in 7-rays 
(all with EGRET). 

2.1.4 The Total Blazar sample 

Combining the three samples yields a total of 126 blazars 
(six of them are present in both the radio and X-ray selected 
samples of BL Lacs), of which 33 detected in 7-rays. We will 
refer to them as the total blazar sample. 

2.2 Multi- frequency Data 

In view of building average SEDs minimizing the bias in- 
troduced by incompleteness, we decided to focus on a few 
well covered frequencies, at which fluxes are available for 
most objects. In a separate paper (Ghisellini et al. 1998) we 
consider a sub-sample of 7-ray loud blazars with extensive 
coverage in frequency with the scope of carrying out detailed 
model fitting for each source. 

We chose the following seven well sampled frequencies, 
that are sufficient to give the basic information on the SED 
shape from the radio to the X-ray band: radio at 5 GHz, 
millimeter at 230 GHz, far infrared (IRAS data) at 60 and 
25 fim, near infrared (K band) at 2.2 /im, optical (V band) 
at 5500 A, and soft X-rays at 1 keV. Data were collected 
from a careful search in the literature and extensive usage of 

the NASA Extragalactic Database (NED)Q. In the radio and 
optical bands the coverage is complete for all the objects in 
the three samples, while unfortunately for mm, far and near 
IR and X-ray fluxes data for some sources are lacking (see 
Table 1). The worse case is the far IR (25 /im) band where 
only 28/126 objects have measured fluxes. 

For each source, at each frequency from radio to optical 
we assigned the average of all the fluxes found in literature. 
Given the large variability these averages were performed 
logarithmically (magnitudes). 

In principle a suitable alternative to the averaging 
would be to consider in each band the maximum detected 
flux (see for instance Dondi & Ghisellini 1995). On one hand 
this choice could be particularly meaningful in view of the 
fact that in the 7-ray band, due to the limited sensitivity 
of detectors, we are biased towards measuring the brightest 
states. On the other hand also this option is biased since 
the value of the maximum flux is strongly dependent on 

t The optical magnitudes have been de-reddened using values of 
Ay derived from the Ab reported in the NED database according 
to the law A v = A B /1.324 (Riecke & Lebofski 1985) 



the observational coverage and for most of the objects we 
only have a few (sometimes a single) observations. More- 
over the strength of this bias is 'band-dependent' and can 
thus significantly affect the determination of the broad band 
spectral shape. As both choices present advantages and dis- 
advantages, and since our goal is a statistical analysis, we 
consider them equally good. The 'averages' option has been 
preferred because it is likely to be more robust with respect 
to the definition of radio-optical SED properties. 

In Table 1 a summary of the collected broad band data 
is reported, with the computed average flux values for each 
object. 

2.2.1 X-ray data 

The knowledge of the X-ray properties is of special rele- 
vance because in this band both the synchrotron and inverse 
Compton processes can contribute to the emission. Since the 
first mechanism is expected to produce a steep continuum 
in this band while the second one should give rise to a flat 
component (a < 1, rising in a vF v plot) the shape of the 
X-ray spectrum can give a fundamental hint for disentan- 
gling the two components and inferring the respective peak 
frequencies. 

We privileged the large and homogeneous ROSAT data. 
base. In fact, a large fraction of the 126 sources (90/126) has 
been observed with the ROSAT PSPC allowing to uniformly 
derive X-ray fluxes and in many cases, that is for 73 targets 
of pointed observations, spectral shapes in the 0.1-2.4 keV 
range (Brunner et al. 1994; Lamer et al. 1996; Perlman et 
al. 1996b; Urry et al. 1996; Comastri et al. 1995, 1997; Sam- 
bruna 1997). X-ray spectral indices were derived from the 
same observation and, when available, we adopted the ax 
resulting from fits with neutral hydrogen column density Nh 
allowed to vary. Some of these 90 objects (17) have been only 
detected in the ROSAT All Sky Survey (RASS) and fluxes 
are published by Brinkmann, Siebert & Boiler (1994) and 
Brinkmann et al. (1995). Monochromatic fluxes (at 1 keV) 
for these sources have been derived from the 0.1 - 2.4 keV 
integrated flux adopting the average spectral index of the 
sample to which they belong (see Table 4) and the value of 
the Galactic column in the source direction (Elvis et al 1989; 
Dickey & Lockman 1990; Lockman & Savage 1995; Murphy 
et al 1996). When more than one observation was available 
we give the average flux. 

Of the remaining 36 sources, 24 belong to the Slew sur- 
vey sample and for them we used directly the Einstein IPC 
flux from Perlman et al. (1996a). The fluxes at 2 keV listed 
by Perlman et al. (1996a) were converted to 1 keV using the 
average ROSAT spectral index of the Slew survey sample 
({ax} =1-40), derived from the 24 sources with a ROSAT 
measured value. 

For other 3 sources, without ROSAT data, we used an 
Einstein IPC flux, bringing the total number of sources with 
measured X-ray flux to 117/126. 

2.2.2 7 -ray data 

Within the three samples only a fraction of blazars were de- 
tected in 7-rays, namely 9/34 in the 1 Jy sample, 8/48 in 
the Slew sample, 20/50 in the FSRQ sample. Four of these 
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Table 2. Basic data for the 30 7-ray detected sources not in- 
cluded in our samples : (1): IAU name; (2): redshift; (3): radio 
flux at 5 GHz; (4): 7-ray flux at 100 McV; (5): EGRET spectral 
index 



(1) 


(2) 


(3) 


( 4 ) 


(5) 


IAU name 


z 


FsGHz 


FiooMcV 








l T \ 

(Jy) 


I T \ 

(nJy) 




0130-171 


1.022 


1.00 


0.122 




0202+149 


1.202 


2.40 


0.383 


1.5 ± 0.1 


0234+285 


1.213 


2.36 


0.296 


1.7 ± 0.3 


0446+112 


1.207 


1.22 


0.470 


0.8 ± 0.3 


0454-234 


1.009 


2.2 


0.143 




0458-020 


2.286 


2.04 


0.364 




0506-612 


1.093 


2.1 


0.062 




0521-365 


0.055 


9.7 


0.139 


1.16 ± 0.36 


0804+499 


1.433 


2.05 


0.322 


1.72 ± 0.38 


0805-077 


1.837 


1.04 


0.404 


1.4 ± 0.6 


0827+243 


0.939 


0.67 


0.226 


1.21 ± 0.47 


0829+046 


0.18 


1.65 


0.132 




0917+449 


2.18 


1.03 


0.075 


0.98 ± 0.25 


1156+295 


0.729 


1.65 


1.727 


1.21 ± 0.52 


1222+216 


0.435 


0.81 


0.278 


1.50 ± 0.21 


1229-021 


1.045 


1.1 


0.250 


1.92 ± 0.44 


1313-333 


1.210 


1.47 


0.098 


0.8 ± 0.3 


1317+520 


1.060 


0.66 


0.079 




1331+170 


2.084 


0.713 


0.091 




1406-076 


1.494 


1.08 


1.013 


1.03 ± 0.12 


1604+159 


0.357 


0.50 


0.260 


0.99 ± 0.50 


1606+106 


1.227 


1.78 


0.312 


1.20 ± 0.30 


1622-297 


0.815 


1.92 


2.416 


1.2 ± 0.1 


1622-253 


0.786 


2.2 


0.336 


1.3 ± 0.2 


1730-130 


0.902 


6.9 


0.258 


1.39 ± 0.27 


1739+522 


1.375 


1.98 


0.236 


1.23 ± 0.38 


1933-400 


0.966 


1.48 


0.158 


1.4 ± 0.2 


2032+107 


0.601 


0.77 


0.192 


1.5 ± 0.3 


2344+514 


0.044 


0.215 


0.8° 




2356+196 


1.066 


0.70 


0.311 





References of Table 2: 

References for the data here reported are listed in the Notes to 
Table 5. 

Notes to Table 2: 

( a > source detected only by Whipple. The given value is the 
integrated flux measured at E > 300 GeV, in units of 10 -11 
photons cm~ 2 s _1 . 

sources (0235+164, 0735+178, 0851+202 and 1652+398) arc 
present in both the BL Lac object samples, giving a net num- 
ber of 7-ray detections of 33 out of 126 blazars. All but one 
of them have been observed by EGRET in the 30 MeV - 30 
GeV band. For 28/32 a 7-ray spectral index has been de- 
termined. One source, 1652+398 (Mkn 501), has only been 
detected at very high energies, beyond 0.3 TeV by ground 
based Cherenkov telescopes (Whipple and HEGRA, Weekes 
et al. 1996; Bradbury et al. 1997), while EGRET yielded 
only an upper limit. It is worth noticing that the detected 
fraction is significantly different between quasars and BL 
Lacs, being respectively 40 ± 10.6 % and 17.1 ± 5.1 % for 
XBLs and RBLs together. However for RBLs only the frac- 
tion detected in 7-rays is 26.5 ± 9.9 %, consistent with that 
of quasars while XBLs only yield 16.7 ± 6.4 %. 

Many other blazars (~ 30) have been detected by 



Table 3. spectral indices used for K-correction of monochro- 
matic fluxes: (1) spectral band; (2) Slew; (3) Uy; (4) FSRQ; (5) 
references. 



band Slew 1 Jy FSRQ refs. 
(1) (2) (3) (4) (5) 



radio 


0.20 


-0.27 


-0.30 


1 


mm 


0.32 


0.32 


0.48 


2 


IRAS 


0.60 


0.80 


1.00 


3 


IR-opt 


0.67 


1.21 


1.52 


4 


X-rays 


1.40 


1.25 


0.83 


5 


7-rays 


0.98 


1.26 


1.21 


5 



References to Table 3: 

(1) Stickel et al. 1994; (2) Gear et al. 1994; (3) derived from IRAS 
data; (4) Bersanelli et al. 1992; (5) this work, see Table 5. 

EGRET, but do not fall in our samples. One can consider 
the group of 7-ray detected objects as a sample in its own 
right, though not a complete one at present, since a sig- 
nificant fraction of the sky has been surveyed though not 
uniformly. This larger sample comprises 66 sources (Fichtel 
et al. 1994; von Montigny et al. 1995; Thompson et al. 1995, 
1997; Mattox et al. 1997 and references therein), of which 60 
have a measured redshift, and 48 an estimate of the spectral 
index. To this set we can add Mkn 501 (already included in 
both our BL Lac samples), and 2344+514, detected only by 
the Whipple telescope (Fegan 1996). We will use this addi- 
tional information to discuss whether the 7-ray properties 
of our samples can be representative of the whole 7-ray loud 
population and if so, to increase the statistics (see section 
3.1). We therefore collected basic data also on all of the 30 
(29 EGRET plus 2344+514) 7-ray detected AGN/blazars 
not included in the complete samples. They are reported in 
Table 2. 

Given the large amount of observations and analysis 
of the same data by different authors, for the selection of 
the flux and spectral index we used the following criteria: 

i) spectral index and flux referring to the same observation, 

ii) data corresponding to a single viewing period, iii) if data 
were analyzed by different authors, the results of the most 
recent analysis are preferred. 

7-ray data are usually given in units of photons cm -2 
s _1 above an energy threshold (e.g. for EGRET E > 100 
MeV). We converted them to monochromatic fluxes at 100 
MeV integrating a power law in photons with the measured 
or assumed (the average) spectral index. 

2.2.3 Luminosities and K-correction 

All fluxes were K-corrected and luminosities were computed 
with the following choices: 

a) we considered lower limits on redshift (4 sources) as de- 
tections, while we assigned the average redshift of the sample 
to the few sources without any estimate (4 in the 1 Jy sam- 
ple, for which (z) = 0.492, and 6 in the Slew survey sample, 
(z) = 0.194); 

b) luminosity distances were calculated adopting Ho=50 
km s _1 Mpc -1 and qo=0; 
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Figure 1. Redshift distributions for the three complete samples 
and the "total blazars sample". Sources detected in the 7-rays 
are indicated by the grey areas. 
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Figure 2. Distributions of L5GHZ for the three complete sam- 
ples and the "total blazars sample". Grey areas indicate 7-ray 
detected objects. 



c) fluxes were K-corrected according to the following pre- 
scriptions. For radio-to-optical data we used average spec- 
tral indices derived from the literature (see Table 3). For 
X-ray and 7-ray data we used measured power law spectral 
indices, when available, or the average index derived for the 
sources of the same sample (see Table 3). 



3 RESULTS 

3.1 Distributions of properties 

Since the fraction of objects detected in 7-rays in the three 
samples is rather small, it is important to ask whether the 
detected sources are representative of each sample as a whole 
or are distinguished in other properties from the rest of the 
objects in it. Moreover we want to verify whether the 7-ray 
detected sources in general differ from those belonging to 
the complete samples. 

We therefore computed the distributions of various 
quantities, i.e. redshift, luminosities and broad band spec- 
tral indices, for objects belonging to each sample. These are 
shown in Figs. 1-5 for the three samples and the total blazar 
one, evidentiating those sources detected in the 7-ray band 
as grey shaded areas in the histograms. The redshift, radio 
(at 5GHz) and X-ray (at IkeV) luminosities (expressed as 
vLv, erg/sec) are shown in Figs. 1-3, while the distributions 
of the broad band spectral indices qro and qrx are plotted 
in Figs. 4 and 5, respectively. 

The redshift distributions (Fig. 1) of the three complete 
samples show the known tendency towards the detection of 
FSRQs at higher z, being the latter ones more powerful radio 
sources, as shown in Fig. 2. 

In the same figure and in Fig. 3, the tendency for XBL 
to have similar X-ray but lower radio luminosities compared 
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Figure 3. Distributions of Lik e y for the three complete sam- 
ples and the "total blazars sample". Grey areas indicate 7-ray 
detected objects. 



to RBL is also apparent. Correspondingly, it appears from 
Figs. 4,5 that oro and «rx increase from XBL to RBL while 
FSRQ have «ro slightly larger and qrx similar to RBLs. 
Later on (section 3.2) we will show that there is a relation- 
ship between these spectral indices and the peak frequency 
of the synchrotron component. 

We note that for all of the distributions there is continu- 
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distinction according to 7— ray detection 



distinction according to inclusion in our samples 





20 
10 



j a 



40 


4— 
41 


42 


43 


44 45 


46 


1 = 
47 


H 














z 1 k^jg^ 














-24 




-23 




-22 








- z 



Figure 4. ctRO distributions for the three complete samples and 
the "total blazars sample". Grey areas indicate 7-ray detected 
objects. 
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Figure 5. orx distributions. As in Fig. 4. 



ity in properties not only between the two BL Lac samples, 
but also between BL Lacs and FSRQs. 

It is clear from Figs. 1-5 that for the RBL and FSRQ 
samples the 7-ray detected sources do not differ from non- 
detected ones in any of the considered quantities while for 
the Slew sample there is a tendency for 7-ray loud sources 
to have larger Lsghz, orx and Oro. This indicates that 
either the radio luminosity or radio-loudness are important 
in determining the 7-ray detection. On the contrary, and in 



Figure 6. Distributions of redshift 2, Lsghzi radio flux Fsghz 
and 7— ray spectral index a 7 for the 7-ray loud sources detected 
by EGRET, where the shaded ones refer only to sources in our 
three complete blazar samples. 



some sense surprisingly, the X-ray luminosity does not seem 
to play an important role with respect to the 7-ray emission, 
although the X-ray band is the closest in energy to the 7- 
rays . We will come back to this issue later on (§3.5). We 
checked the possible difference of means and variances of 
the distributions with the t-student's test and only for the 
cvrx of the Slew sample the significance is higher than 95 
per cent. 

Except for the case of the Slew survey, we conclude that 
the 7-ray detected sources are representative of the samples 
as a whole, being indistinguishable from the others in terms 
of radio to X-ray broad band properties and power. 

We also checked that the 7-ray detected sources belong- 
ing to our samples are homogeneous with respect to all of the 
7-ray blazars detected so far. In Fig. 6 we compare the red- 
shift distributions, the radio luminosities and fluxes and the 
7-ray spectral index. The grey shaded areas represent the 
7-ray sources belonging to the complete samples considered 
here. We conclude that there is no significant difference. 

Nevertheless, we are aware that the limited sensitivity 
of the EGRET instrument implies that at a given radio flux, 
only the 7-ray loudest sources are detected. Therefore the 
non detected ones are probably on average weaker in 7-rays. 
Impey (1996) quantified this effect by taking into account 
the correlation between radio and 7-ray luminosities (see 
§3.3.1), and other observables. Assuming a Gaussian dis- 
tribution of the 7-ray to radio flux ratio he estimated the 
width of the distribution and the "true" ratio referring to 
the whole population, which could be a factor 10 lower than 
the observed one. There could be a real spread in the intrin- 
sic properties of the blazar population, the 7-ray detected 
blazars being intrinsically louder than the rest of the pop- 
ulation. Alternatively this may be due to variability, i.e. a 
source is detected only when it undergoes a flare. The ob- 
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served ratio would thus refer to flaring states, while the av- 
erage level of each source would be lower. 

Since variability is a distinctive property of blazars and 
has been observed to occur also in 7-rays, often with ex- 
tremely large amplitude (greater than a factor 10) (e.g. 3C 
279, Maraschi et al 1994; PKS 0528+134, Mukherjee et al. 
1996; 1622-297, Mattox et al. 1997) the latter alternative is 
likely although the problem remains an open one. We con- 
clude that the average 7-ray luminosities computed here 
are necessarily overestimated. However we chose not to cor- 
rect for this effect given the uncertainties. In particular the 
"bias factor" for different classes of blazars could be differ- 
ent if their 7-ray variability properties (amplitude and duty 
cycle) are (e.g. Ulrich, Maraschi & Urry 1997). 



3.2 Synchrotron peak frequency 

As previously noted the SEDs clearly show a broad peak 
between radio and UV-X-rays. In order to determine the 
position of the peak of the synchrotron component in in- 
dividual objects with an objective procedure, we fitted the 
data points for each source (in a v vs. uL„ diagram) with a 
third degree polynomial, which yields a complex SED pro- 
file, with an upturn allowing for X-ray data points not to 
lay on the direct extrapolation from the lower energy spec- 
trum. In many cases there is evidence that the X-ray com- 
ponent, even in the soft ROSAT PSPC band, is due to the 
inverse Compton process (e.g. Sambruna 1997; Comastri et 
al. 1997). Thus to impose that the X-ray point smoothly 
connects to the lower energy data, as would happen in a 
parabolic fit, could be misleading for a determination of 
the synchrotron peak frequency. We used a simple parabola 
when the cubic fit was not able to find a maximum, which 
typically happens when the peak occurs at energies higher 
than X-rays. In fact when the peak moves to high enough 
frequencies (typically beyond the IR band), the X-ray flux 
is completely dominated by the synchrotron emission, and 
the results given by the cubic and parabolic fits are fully 
consistent. In 8 cases neither the cubic nor the parabolic fit 
were able to determine a peak frequency /luminosity mainly 
due to the paucity of data points. 



3.2.1 Synchrotron Peak Frequency vs. Luminosity 

The peak frequencies derived with the above procedure (de- 
fined as the frequencies of the maximum in the fitted poly- 
nomial function) are plotted in Fig. 7a,b,c versus the radio 
and 7-ray luminosities and vs the corresponding peak lu- 
minosities, as determined from the fits. Let us stress once 
again the continuity between the different samples. Consid- 
ering the samples together strong correlations are present 
between these quantities, in the sense of f pC ak,sync decreasing 
with increasing luminosity. The results of Kendall's r sta- 
tistical test (Table 4) show that the correlations are highly 
significant. 

Since on one hand in flux limited samples spurious cor- 
relations can be introduced by the luminosity /redshift rela- 
tion and on the other hand the correlations might be due to 
evolutionary effects genuinely related to redshift, we checked 
its role in two ways. We estimated the possible correlation of 
the relevant quantities with redshift directly, and performed 



48 



46 



--4 



Of 

o 



42 - 



40 - 



50 



4ii 



O 



44 



48 - 



46 



O 



42 



40 - 



@5 GHz 




12 14 


16 18 


log 


(v , ) 
^ peak, sync/ 


_ 1 1 1 1 1 1 
A 


1 1 1 1 1 1 1 1 1 1 

@y — rays 




A 






Q A 




^A A 




A 


® 






O 


® 




O 






: (b) 

~ , 1 , , , 1 


, 1 , , , 1 , , , ~ 


12 14 


16 18 


log 


iy u ) 

^ peak, sync' 



peak, sync 




(c) 



A FSRQ 

O BLL (Uy) 

® BLL (Slew) 



12 14 16 

1op(i> , ) 

to V peak, sync-' 



1 s 



Figure 7. The peak frequency of the synchrotron component, 
v peak,sync! as derived with the polynomial fits, plotted against 
a) the radio luminosity L5GHZ1 b) the 7-ray luminosity L 7 , and 
c) the fitted peak luminosity of the synchrotron component, 

Lpeak,sync- 



8 G. Fossati, L. Maraschi, A. Celotti, A. Comastri, G. Ghisellini 



0.8 



0.6 



0.4 



0.3 




■ radio-X 
□ radio-optical 



□ 



1 2 



14 16 

\og{v 



peak, sync/ 



113 Hz, with a RO = 0.76 and a RX = 0.85, 
10 16 Hz (a RO = 0.35 and a RX = 0.56). 



Figure 8. The broad band spectral indices a R o and a R x are 
plotted vs t'peak^ync- The curved lines overlayed to the data 
points represent the relations a R o-f pe ak,syiic and "RX-fpcak.sync 
obtained from a "synthetic" set of SEDs. Details on the adopted 
analytical parameterization are given in text, §4. Two examples 
of typical SEDs are reported for reference in the inset: the top 
one peaking at v ~ 10 1 
the bottom one at v 
The three vertical lines mark the frequencies corresponding to 5 
GHz, 5500 A, and 1 keV, entering in the definitions of a R o and 
o R x- 



partial correlation tests between two quantities subtracting 
out the common dependence on z (Padovani 1992b) (see re- 
sults in Table 4). In addition, in order to have an indepen- 
dent check on the redshift bias, we also considered the sig- 
nificance of the correlations restricting them to objects with 
z < 0.5 (see Table 4). 

The correlation between ^ pea k,sync and Lsghz still holds 
after subtraction of the very strong dependence on redshift. 
The same is true for the relation between f pC ak,sync and the 
7-ray luminosity, although the significance is much smaller, 
due to the smaller number of sources. On the other hand 
the correlation between i/peak.sync and L pea k,sync is strongly 
weakened when subtracting the redshift effect. 

Considering only the z < 0.5 interval the significance of 
the first two correlations persists and does not change when 
the redshift dependence is subtracted. These values can then 
be considered as irreducible, being the signature of a true 
dependence of ^ pe ak,sync on luminosity. This result can also 
be read as a check of the reliability of the partial correla- 
tion procedure. On the contrary the correlation i^ pca k,sync vs. 
L P cak,sync disappears at low redshifts, due to the narrow in- 
terval of values spanned by L poak>sync , that varies less with 
the change of peak frequency than both radio and 7-ray 
luminosity do. 
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3.2.2 Synchrotron Peak Frequency vs. Broad Band 
Spectral indices 

The relations between the synchrotron peak frequency and 
each of the two point spectral indices a R ,o and a R x are 
shown in Fig. 8. Also these quantities are strongly correlated 
(see Table 4). In fact recent papers (e.g. Maraschi et al. 1995; 
Comastri et al. 1995; Comastri et al. 1997), suggested that 
the position of the synchrotron peak could be devised from 
the values of broad band spectral indices. 

We see that the knowledge of any of the two spec- 
tral indices is enough to guess the position of the peak 
of the synchrotron component, except for some ranges, 
namely f pC ak,sync > 

and i/pcak.sync < 10 

explained bearing in mind the typical shape of the blazar 
SEDs (see inset in Fig. 8): when the spectrum peaks at low 
frequencies, X-rays are typically dominated by the inverse 
Compton, flat spectrum, component whose luminosity level 
is strongly correlated with the radio one (Fossati et al. 1997), 
and then the X-ray/radio ratio (i.e. a R x) tends to a fixed 
value. Conversely the Compton component begins to dom- 
inate the (ROSAT) X-ray band when a R x ~ 0.75, corre- 
sponding to fp C ak,sync < 3 x 10 14 Hz. It is interesting to note 
that the adopted dividing threshold between LBL and HBL 
has been set to this same value from purely practical pur- 
poses, while in the light of the result above it assumes a more 
"physical" meaning. LBL sources have Compton-dominated 
soft-X-ray emission, while in HBL this is pure synchrotron. 

At the other end of the spectrum a problem arises when 
t'peak.sync moves at energies higher than that used to com- 
pute the broad band spectral index. The reason is that the 
ratio between, for instance, optical and radio luminosity 
is no longer sensitive to the peak moving further towards 
higher frequencies, because both the radio and optical bands 
lay on the same (rising) branch of the synchrotron "bump" . 

For comparison we draw in Fig. 8 the loci of a R o- 
fpcak.sync and aRx-v peak, sync obtained from a set of SEDs 
of the kind reported in the inset, and that we are going to 
discuss in more detail in section 3.5. The parameterization 
describes the observed features very well. 



3.2.3 Synchrotron Peak Frequency vs. "/-ray dominance 

In Fig. 9 ,sync is plotted against the 7-ray dominance 
parameter, defined as the ratio between the 7-ray and the 
synchrotron peak luminosities. A strong correlation (see Ta- 
ble 4) is present over four orders of magnitude in i7>cak,sync, 
in the sense of a decrease in the 7-ray dominance with an 
increase of the synchrotron peak frequency. In the same fig- 
ure we plotted also the ratio between the 7-ray and optical 
luminosities, to check if the latter could eventually be a good 
indicator of the 7-ray dominance, with the advantage of be- 
ing an observed quantity. In fact there is little difference, 
at most a factor 3 for a quantity spanning more than three 
decades. 



3.3 Average SEDS 

Having discussed extensively the possible biases introduced 
by the limited number of 7-ray detected sources in the com- 
plete samples we construct here the average SEDs for each 
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Table 4. Correlation probabilities according to the Kendall's r test. (1-2): quantities considered in the correlation; (3-10): significance of 
the "null-hypothesis", i.e. that the correlation is the result of pure chance, for the various quantities: in columns (3-6) for tests performed 
without any redshift constraint, while in (7-10) results for tests taking into account only sources with z > 0.5. In (3), (7) the "face value" 
correlation between xi and X2. In (4—5) and (8—9) the significance of the correlation of each quantity xi and X2 with redshift. In (6) 
and (10) the "net" xi/x2 correlation remaining from (3) and (7) after subtraction of the redshift dependence of xi and X2 via partial 
correlation algorithm. A dash is reported the correlation is not significant (i.e. probability < 90 %. 

z unconstrained (# 109°) z < 0.5 (# 51°) 

XI X2 Xl/X2 xi /z X2/z xi/x2—Z Xl/X2 xi/z X2/z xi/x2— Z 

(1) (2) (3) (4) (5) (6) J7) (8) (9) (10) 

1.3e-8 [ — ] [4.0c-4] 9.2e-9 

[ - ] [1.2e-7] — 

1.8c-7 [ — ] [ — ] 2.9o-7 

4.1c-15 [ — ] [ — ] 8.4c-15 



Correlations involving 7-ray data 



EGRET sources included in our samples 

z unconstrained (# 31) z < 0.5 (# 11) 



^peak^ync ^-"y 


2.1c-3 


[3.4e-2] 


[1.5e-10] 


1.8e-2 


2.4e-2 


[ - ] [2-4e-3] 


6.4e-2 


LsGHz L 7 


6.0c-ll 


[5.5e-10] 


[1.5e-10] 


4.2e-5 


3.9e-3 


[5.2e-2] [2.4e-3] 


2.4e-2 


^peak^sync /Lpcak.sync 


2.7c-5 


[3.4c-2] 


[1.5e-5] 


2.2e-4 


1.4e-3 


[ - ] [7.36-1] 


4.1c-3 


fpcak.sync L t/ L 5500A 


2.5e-6 


[3.4e-2] 


[4.2e-5] 


1.7e-5 


2.4e-3 


[-] [-] 


6.8c-3 



sources in the whole EGRET sample 

z unconstrained (# 60) z < 0.5 (# 15) 



LsGHz L 7 4.8e-15 [4.8e-15] [1.8e-14] 2.4e-6 1.8e-3 [l.le-2] [1.4c-4] 4.1e-2 



"peak.sync LsGHz l-2c-16 [l.lc-9] [1.9e-31] l.le-11 

lpcak.sync Lpeak.sync 4.9C-8 [l.lc-9] [9.5C-27] 6.5e-2 

^pcak, S ync o R o 11^23 [l.lc-9] [7.7c-ll] 4.7e-19 

lpcak.sync «rx 4.2c-17 [l.lc-9] [6.0c-5] 1.9c-14 



Notes to Table 4: 

(°) wc considered only sources with at least a lower limit on redshift, and for which it has been possible to determine the "synchrotron" 
peak frequency by means of the polynomial fit. 



sample. We will come back later to the bias problem (Section 
4). 

The averaging procedure has been performed on the 
logarithms of the luminosities at each frequency. Apart from 
the problems in the 7-ray range discussed above the incom- 
pleteness of the data coverage at some frequencies could also 
introduce a bias in the average values. For instance in the 
Slew survey sample only 10/48 objects have a flux measured 
at 230 GHz, and they are the more luminous sources at 5 
GHz. Averaging independently L230GHZ (for 10 objects) and 
LsGHz (for 48 objects) we would obtain a ratio between the 
two luminosities higher than that derived considering only 
the subsample of 10 sources, and presumably higher than 
the actual one, too. 

To reduce this kind of effect we first normalized the 
monochromatic luminosities to the radio luminosity for each 
source, we computed average ratios (log(Z/„* /I/5GHz))|sub, 
considering only the subsample of sources with a measured 
flux at v* , and used that ratio to compute the average 
monochromatic luminosity at v* for all sources in the sample 



as (log(L„»))|aii = (lOg(L 5 GHz))|aii + (log^. /L SG Hz )) \sub ■ 

In this way we basically averaged the spectral shape between 
is* and 5 GHz for the measured objects and assigned that 
spectral shape to the sample. 

The X-ray and 7-ray spectral indices have been aver- 
aged with a simple mean, without weighting. 

The average broad band spectra for each of the three 
samples are shown in Fig. 10. The 6 sources common to the 
radio and the X-ray selected BL Lac samples are considered 
in both of them. Average luminosities entering Fig. 10 arc 
reported in Table 5 together with the number of sources 
contributing at each frequency. 

It is apparent from Fig. 10 that the three samples re- 
fer to objects with different average integrated luminosities 
and that the peak frequency of the power emitted between 
the radio and the X-ray band moves from the X-ray to the 
far infrared band going from the XBL to the FSRQ sam- 
ples as anticipated from the analysis of single objects in the 
previous section 3.2. Correspondingly the 7-ray luminosi- 
ties increase and the 7-ray spectra steepen suggesting that 
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Figure 9. The 7-ray dominance (according to two definitions, 
see text) versus the synchrotron peak frequency fpeaj^sync 
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Figure 11. X-ray and 7-ray spectral indices plotted against ra- 
dio luminosity. 
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Figure 10. The average SEDs for each of the samples are shown. 
From top to bottom (referring to radio luminosity) Wall & Pea- 
cock FSRQs (empty boxes), 1 Jy BL Lac sample (filled boxes) 
and Slew survey BL Lac sample (triangles) . These latter two are 
in reversed order in X-ray band, the lowest spectrum being that 
of the 1 Jy sample. 



also the peak frequency of the high energy emission moves 
to lower frequencies. The overall similarity and regularity of 
the SEDs of the different samples as well as the continuity 
in the properties of the individual objects discussed in sec- 
tion 3.2 suggest a basic similarity of all blazars irrespective 
of their original classification and different appearance in a 
specific spectral band. 

We therefore considered the merged total sample with 
the scope of finding the key parameter(s) governing the 
whole blazar phenomenology. Since luminosity appears to 



have an important role in that it correlates with the main 
spectral parameters we decided to bin the total blazar sam- 
ple according to luminosity, irrespective of the original clas- 
sification. We used the 5 GHz radio luminosity which is 
available for all objects. It may be desirable to use the to- 
tal integrated luminosity which in all cases is close to the 
7-ray one. However the latter is only available for a few ob- 
jects. We note that a correlation between 7-ray and radio 
luminosity has been claimed by many authors using different 
techniques (Dondi & Ghisellini, 1995; Mattox et al. 1997). It 
is however still being debated whether it is true or it arises 
from selection effects, connected with the common redshift 
dependence of luminosities, and with the exclusion of up- 
per limits, which could favour the appearance of a spurious 
correlation. It is worth mentioning that Mucke et al. (1997) 
using a technique designed to take into account both these 
effects did not find any significant correlation between ra- 
dio and 7-ray data for a sample of 38 extragalactic EGRET 
sources. 

We also checked this correlation on both the 31 EGRET 
detected sources included in our samples and the larger 
"comparison sample" of 62 EGRET sources. In Table 4 
we report the significance of the correlation, together with 
its value after subtracting the common redshift dependence 
through a partial correlation test, and its significance for 
samples restricted to z < 0.5. In all cases the radio and 
7-ray luminosities correlate significantly. 

In Figs. 11 Qx and a 7 for individual sources are plotted 
against the radio power, both showing a good correlation 
with it. Comastri et al. (1997) discussed the interesting con- 
sequences of the apparent anti-correlation between X-ray 
and 7-ray spectral indices, without relating it to any "abso- 
lute" parameter, such as luminosity. Here again we see that 
these other spectral properties have a dependence on radio 
luminosity. 

Since in some luminosity bins the number of 7-ray de- 
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Table 5. Average luminosities, for each sample and for the total one (divided in bins of radio luminosity) 
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Complete Samples Total Blazar Sample (log(L 5 Gjjz) intervals) 



Band Slew 1 Jy W&P <42 42-43 43-44 44-45 >45 



5 GHz 41.71 43.69 

48 34 

230 GHz 43.11 45.01 

10 34 

60 iiui 44.17 45.94 

12 19 

25 inn 44.25 46.07 

10 15 

K-band 44.64 45.86 

23 31 

V-band 44.91 45.68 

48 34 

1 keV 44.94 44.72 

48 32 

100 MeV 44.45 46.50 

7 8 

a x 1-40 ± 0.07 1.25 ± 0.09 

24 31 

a 7 0.98 ± 0.32 1.26 ± 0.26 

6 6 



44.81 41.24 42.47 

50 38 10 

46.11 42.64 43.77 

50 5 7 

46.84 43.73 44.65 

13 6 5 

46.93 43.74 44.95 

8 4 6 

46.49 44.42 45.04 

28 13 10 

46.58 44.61 45.01 

50 38 10 

45.98 44.81 44.11 

43 38 10 

47.93 44.24 44.79 

20 3 5 

0.83 ± 0.08 1.37 ± 0.09 1.55 ± 0.15 

24 16 8 

1.21 ± 0.09 0.64 ± 0.07 0.73 ± 0.47 

18 2 4 



43.71 44.54 45.39 

17 44 17 

45.13 45.83 46.63 

15 44 17 

46.09 46.65 47.61 

8 16 2 

46.08 46.79 47.69 

7 9 2 

45.96 46.27 47.21 

15 32 6 

45.82 46.27 47.21 

17 44 17 

44.92 45.66 46.50 

15 42 12 

46.67 47.71 48.68 

9 33 12 

1.16± 0.14 1.11 ± 0.08 0.57 ± 0.13 

14 26 9 

1.37 ± 0.31 1.06 ± 0.13 1.30 ± 0.08 

7 25 11 



References to Table 5: 

(5 and 230 GHz): Becker, White & Edwards 1991; Bloom ct al. 1994; Gear ct al. 1986; Gear 1993a; Gear ct al. 1994; Kiihr et al. 1981; 
Kiihr & Schmidt 1990; Perlman ct al. 1996a; Reuter ct al. 1997; Steppe ct al. 1988, 1992, 1993; Stevens et al. 1994; Stickel et al. 1991; 
Stickel ct al. 1993; Stickel, Meisenheimer & Kiihr 1994;Tornikovski et al. 1993, 1996; Terasranta ct al. 1992. 

(IR-optical data): Allen, Ward & Hyland 1982; Ballard et al. 1990; Bersanelli et al. 1992; Bloom ct al. 1994; Brindle et al. 1986; Brown 
et al. 1989; Elvis et al. 1994; Falomo ct al. 1988; Falomo et al. 1993a; Falomo ct al. 1993b; Falomo, Scarpa & Bersanelli 1994; Gear ct 
al. 1986; Gear 1993b; Glass 1979, 1981; Holmes et al. 1984; Impcy & Brand 1981; Impcy & Brand 1982; Impey et al. 1982; Impey et al. 
1984; Impey & Neugebauer 1988; Impey & Tapia 1988, 1990; Jannuzi, Smith & Elston 1993, 1994; Landau ct al. 1986; Lepine, Braz & 
Epchtein 1985; Lichtfield et al. 1994; Lorenzetti et al. 1990; Mead et al. 1990; O'Dell et al. 1978; Pian ct al. 1994; Sitko & Sitko 1991; 
Smith et al. 1987; Stevens et al. 1994; Wright, Abies & Allen 1983. 

(X— rays): Brinkmann et al. 1994; Brinkmann et al. 1995; Brunner et al. 1994; Comastri et al 1995; Comastri et al 1997; Lamer, 
Brunner & Staubert 1996; Maraschi ct al. 1995; Perlman ct al. 1996a; Perlman et al. 1996b; Sambruna 1997; Urry et al. 1996. 
(7-rays): Bertsch et al. 1993; Catanese et al. 1997; Chiang et al. 1995; Dingus ct al. 1996; Fichtel et al. 1994; Hartman ct al. 1993; Lin 
et al. 1995; Lin ct al. 1996; Madejski et al. 1996; Mattox ct al. 1997; Mukherjee et al. 1995, 1996; Nolan ct al. 1996; Quinn ct al. 1996; 
Radecke et al. 1995; Shrader ct al. 1996; Srcckumar ct al. 1996; Thompson et al. 1993; Thompson et al. 1995; Thompson et al. 1996; 
Vestrand, Stacy & Srcckumar 1995; von Montigny et al. 1995. 



tected sources is small, we used an indirect procedure to 
associate 7-ray fluxes and spectra to our average SEDs, tak- 
ing advantage of the whole body of information regarding 
the 7-ray properties of blazars. Namely for each luminosity 
bin (L 7 ) and (a 7 ) were computed from blazars in the gen- 
eral EGRET-detected sample falling into the same Lsghz 
bin. The basic assumption is the uniformity of the spectral 
properties, as discussed in section 3.1. The resulting SEDs 
are shown in Fig. 12 and average luminosities, X-ray and 
7-ray spectral indices, and number of sources are reported 
in Table 5. The most interesting result is that the trends 



pointed out for the three separate sub-classes of blazars 
(Fig. 10) hold for the total blazar sample, irrespective of 
the original classification of sources, when the radio lumi- 
nosity is adopted as the key parameter characterizing each 
object. 



4 DISCUSSION 

In Fig. 12 we superimposed to the averaged data a set of 
(dashed) lines, whose main goal is to guide the eye. The 
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Figure 12. Average SEDs for the "total blazar sample" binned 
according to radio luminosity irrespective of the original classifica- 
tion. The ovcrlayed curves are analytic approximations obtained 
according to the one-parameter-family definition described in the 
text. 

radio-X-ray SED is approximated with a power law starting 
in the radio domain continuously connecting at v ~ 5 x 10 11 
Hz with a parabolic branch. This latter describes the peak of 
the SED and its steepening beyond it. In soft X-rays a rising 
power law, representing the onset of the hard inverse Comp- 
ton component is summed to this curved "synchrotron" 
component. The normalization of this second X-ray compo- 
nent is kept fixed relative to the radio one. Based on our find- 
ings (see Fig. 7a), we then assume that the peak frequency 
of the synchrotron spectral component is (inversely) related 
to radio luminosity. The simplest hypothesis of a straight 
unique relationship between t'pcak.sync and Lsghz does not 
give a good result when compared with the average SEDs. 
We then allow for a different SED-shape/luminosity depen- 
dence for high and low luminosity objects, a distinction that 
turns out to roughly correspond also to that between objects 
with and without prominent emission lines. We adopted a 
"two-branch" relationship between ^ pC ak,sync and L 5 ghz in 
the form of two power laws ^ poa k,sync oc L 5( ? Hz , with 77 = 0.6 
or n — 1.8 for log(LsGHz) higher or smaller than 42.5, re- 
spectively. The shape of the analytic SEDs is parabolic with 
a smooth connection to a fixed power law in the radio and 
the loci of the maxima as defined above. A full description 
of the parameterization can be found in Fossati et al. (1997) 
where a similar scheme was proposed to account for the 
source number densities of BL Lacs with different spectral 
properties (LBL and HBL). 

The analytic representation of the second spectral com- 
ponent (X-ray to 7-rays) is a parabola of the same width as 
the synchrotron one, and has been obtained assuming that: 
(a) the ratio of the frequencies of the high and low energy 
peaks is constant (f pC ak,Com P t/f peak, sync — 5 x 10 8 ), (b) the 
high energy (7-ray) peak and radio luminosities have a fixed 
ratio, i^Lpcak.gamma/V'iGHzLsGHz ~ 3 x 10 3 . Given the ex- 
treme simplicity of the latter assumptions, it is remarkable 



that the phenomenological analytic model describes the run 
of the average SEDs reasonably well. The worst case refers 
to the second luminosity bin: the analytic model predicts a 
7-ray luminosity larger than the computed bin average by 
a factor of 10 (but predicts the correct spectral shape). We 
note that only 5 7-ray detected objects fall in this bin. 

The results derived from the above analysis (see in par- 
ticular Figs. 10-12) can then be summarized as follows: 

(i) two peaks are present in all the SEDs. The first one 
(synchrotron) is anticorrelated with the source luminosity 
(see Figs. 7 and Table 4), moving from ~ 10 16 - 10 17 Hz 
for less luminous sources to ~ 10 13 — 10 14 Hz for the most 
luminous ones. 

(ii) the X-ray spectrum becomes harder while the "y-ray 
spectrum softens with increasing luminosity, indicating that 
the second (Compton) peak of the SEDs also moves to lower 
frequencies from ~ 10 24 — 10 25 Hz for less luminous sources 
to ~ 10 21 — 10 22 Hz for the most luminous ones; 

(iii) therefore the frequencies of the two peaks are corre- 
lated: the smaller the i^ pC ak.sync the smaller the peak fre- 
quency of the high energy component; a comparison with 
the analytic curves shows that the data are consistent with 
a constant ratio between the two peak frequencies; 

(iv) increasing Lsghz increases the 7-ray dominance, i.e. 
the ratio of the power emitted in the inverse Compton and 
synchrotron components, estimated with the ratio of their 
respective peak luminosities (see also Fig. 9). 

The fact that the trends present when comparing the 
different samples (e.g. Fig. 10), persist when the total blazar 
sample is considered and binned according to radio lumi- 
nosity only, suggests that we deal with a continuous spectral 
sequence within the blazar family, rather than with sepa- 
rate spectral classes. In particular the "continuity" clearly 
applies also to the HBL - LBL subgroups: HBL have the 
lowest luminosities and the highest peak frequencies. 

An interesting result apparent from the average SEDs is 
the variety and complexity of behaviour shown in the X-ray 
band. As expected the crossing between the synchrotron and 
inverse Compton components can occur below or above the 
X-ray band affecting the relation between the X-ray lumi- 
nosity and that in other bands. A source can be brighter than 
another at 1 keV being dimmer in the rest of the radios- 
ray spectrum except probably in the TeV range. This effect 
narrows the range of values spanned by Li kc v and explains 
why 7-ray detected sources do not select a particular range 
in the X-ray luminosity distributions (see Fig. 3) while this 
happens for Lsghz- 

Using this simple scheme of SED parameterization we 
can compute the luminosities in the EGRET (30 MeV - 3 
GeV) and Whipple (0.3 - 10 TeV) bands. These are plotted 
in Fig. 13 (bottom panel) together with their ratio with the 
radio and X-ray luminosities (top and middle panel respec- 
tively). 

It is easy to recognize that for a given radio flux sources 
with t'peak.sync around 10 14 have the largest relative flux in 
the EGRET band, because the peak of the Compton compo- 
nent falls right there (Fig. 13, top panel). For higher values 
of i-'peak.sync the 7-ray peak moves to higher energies too 
and the contribution in the EGRET band is reduced. For 
sufficiently high ^ pe ak, S ync the 7-ray peak reaches the TeV 
band where it becomes detectable. 
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Figure 13. The predicted relations with ^ pea k,sync of the follow- 
ing quantities: the ratio of the EGRET (dashed line) and Whipple 
(solid line) fluxes (luminosities) with the radio one (top panel); 
the ratio of the EGRET and Whipple fluxes (luminosities) with 
the X— ray one (middle panel); EGRET and Whipple absolute lu- 
minosities (bottom panel). The 7-ray luminosities are integrated 
in a band approximately corresponding to that of EGRET or 
Whipple telescopes, 30 MeV-3 GeV and 0.3-10 TeV, respectively. 

Qualitatively the same general behaviour is present also 
in the ratios between EGRET/ Whipple fluxes and the X- 
ray one (Fig. 13, middle panel). There are however a cou- 
ple of significant differences: firstly the EGRET/X-ray ratio 
profile, while still peaking around 10 14 Hz, is sharper than 
in the EGRET/radio case, secondly the TeV relative flux 
distribution is broader and skewed towards lower values of 
the synchrotron peak frequency. Thus, for a given X-ray 
flux (as would be the case in a flux limited X-ray selected 
sample) only those sources falling in the restricted interval 
fpeak.sync ~ 10 13 - K) 15 Hz would have a flux in the EGRET 
band high enough to be detectable. On the other hand, for 
the TeV band it turns out that the chance of being observ- 
able is not confined to very extreme HBLs, with X-ray syn- 
chrotron peaks, but also intermediate BL Lacs could reach 
a comparable TeV flux. 

Since ^peak.sync is directly related to both oro and qrx 
we can understand now the tendency (section 3.1) of 7-ray 
detected sources in the Slew survey to have larger values of 
arto an d chrx- Moreover, due to the fact that in the Slew 
sample LBLs are only a small fraction, the discussion above 
explains also the lower EGRET detection rate with respect 
to other blazar samples. 

The proposed scenario relates the shape of the contin- 
uum to the total source power. It follows that predictions 
of this unifying scheme on both the detectability of blazars 
at 7-ray energies (in view of more sensitive 7-ray detectors, 
e.g. GLAST, improved Cherenkov telescopes, etc.), and their 
contribution to the 7-ray diffuse background depend on the 
combined effects of the SED shape, the luminosity functions 
and possibly evolution (Fossati et al., in preparation; see 
also Stecker, de Jager & Salamon 1996). In particular, an 



interesting and testable prediction of the scheme is the ab- 
sence of high luminosity sources with synchrotron peaks in 
the X-ray range and strong associated TeV emission. 



4.1 Interpretation 

The extreme "regularity" of the SEDs of blazars and in 
particular the trends discussed above must derive from the 
common underlying physical processes. The common sce- 
nario envisaged is that of a relativistic jet pointing close 
to the line of sight. Assuming the simple case of a single 
(homogeneous) zone model the shape of the SED depends 
on the spectrum of the high energy electrons radiating via 
synchrotron and inverse Compton, the magnetic field and 
the nature of seed photons for the inverse Compton pro- 
cess. The latter could be the synchrotron photons themselves 
(synchrotron self Compton, SSC) or photons outside the jet 
("external Compton" EC). In the following we discuss qual- 
itatively the implications of the suggested trends for the two 
scenarios. 

Let us assume a constant bulk Lorentz factors in all 
blazars. Should the (homogeneous) SSC model be valid for 
all sources, it is easy to see that the (approximately) con- 
stant ratio between the high and low peak frequencies yields 
an (approximately) constant value for the energy of the par- 
ticles radiating at the peaks (e.g. Ghisellini, Maraschi, Dondi 
1996). If the energy of the radiating particles is similar in all 
sources the different peak frequencies should result from a 
systematic variation in magnetic field strength, HBLs hav- 
ing the highest, FSRQs the lowest random field intensity. 

Should instead the soft photons upscattered to the 7- 
ray range be produced outside the jet at a "typical" fre- 
quency (the same for all sources) the condition of a constant 
ratio between the peak frequencies implies a constant value 
of the magnetic field (Ghisellini, Maraschi, Dondi 1996). As 
a consequence the energy of the particles radiating at the 
peaks should vary along the spectral sequence being lower 
in FSRQs and higher in HBLs. 

It could also be that there is a smooth transition be- 
tween the SSC and EC mechanisms depending on the phys- 
ical conditions outside the jet. In all cases however the role 
of the luminosity, which phenomenologically appears domi- 
nant, does not find an immediate physical justification, al- 
though one could find plausible arguments to link it to the 
parameters mentioned above and in particular to the condi- 
tions surrounding the jet. 

In a separate paper (Ghisellini et al. 1998), we perform 
model fits to the spectral energy distributions of 51 indi- 
vidual objects, deriving the (model dependent) physical pa- 
rameters for each source. These computations indeed suggest 
the idea that the blazar sequence follows from a transition 
from the SSC to the EC scenario, RBLs being the inter- 
mediate objects. The computed radiation energy densities, 
which determine the amount of radiative cooling, increase 
with increasing source luminosity and may be responsible 
for the lower energy of the particles radiating at the peaks 
in higher luminosity sources. 

The likely possibility that the external photon field in- 
volved in the EC process is (or is related to) the radiation re- 
processed as broad emission lines, seems to be at least qual- 
itatively in agreement with the observational evidence con- 
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cerning the emission line luminosity in the suggested blazar 
sequence. 



5 CONCLUSIONS 

The main conclusion of this work is that despite the dif- 
ferences in the continuum shapes of different sub-classes of 
blazars, a unitary scheme is possible, whereby Blazar con- 
tinua can be described by a family of analytic curves with 
the source luminosity as the fundamental parameter. The 
"scheme" (admittedly empirical) determines both the fre- 
quency and luminosity of the peaks in the synchrotron and 
inverse Compton power distributions (and therefore also the 
7-ray luminosity in the EGRET range) starting from the 
radio luminosity only. The main suggested trend is that 
with increasing luminosity both the synchrotron peak and 
the inverse Compton peak move to lower frequencies and 
that the latter becomes energetically more dominant. The 
scheme is testable, for instance it predicts that sources emit- 
ting strongly in the TeV band have relatively low intrinsic 
luminosity. 

The "spectral sequence" finds a plausible interpretation 
in the framework of relativistic jet models radiating via the 
synchrotron and inverse Compton processes if the physical 
parameters (magnetic field and/or critical energy of the ra- 
diating electrons) vary with luminosity or if photons outside 
the jet become increasingly important as seed photons for 
the inverse Compton process in sources of larger luminosity. 
The latter alternative is supported at least qualitatively by 
the increasing dominance of emission lines in higher lumi- 
nosity objects. 

The proposed scenario, in which the intrinsic jet power 
regulates, in a continuous sequence, the observational prop- 
erties from the weaker HBL, through LBL, to the most pow- 
erful FSRQs, also fits in very nicely with the unification of 
FR I and FR II type radio galaxies as proposed by Bicknell 
(1995). After a long debate the prevailing view is that FR I 
and FR II radio galaxies both contain relativistic jets which 
can be decelerated giving rise to the FR I morphology de- 
pending on the kinetic power in the jet and the pressure of 
the ambient medium. 

The whole radio-loud AGN population could be unified 
in a two parameter space one being the intrinsic jet power, 
the other the viewing angle. An interesting point for future 
discussion is whether a third parameter associated with the 
luminosity of an accretion disk is necessary or is already 
implicitly and uniquely linked to the jet power. 
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wp 
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0.286 
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0.883 










0145+138 


Slew 


0.125 
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0.527 


0.267 




0.162 


0.90 ± 0.40 


0502+675 


Slew 




0.028 












0.744 


3.903 








0507-040 


Slew 


0.304 


0.022 












0.101 


2.306 








0528+134 


wp 


2.060 


0.930 


1 


308 








0.375 


0.951 


0.54 ± 0.29 


1.738 


1.60 ± 0.10 


0537-441 


Uy 


0.896 


2.110 


2 


291 


0.495 


0.215 


8.841 


1.626 


0.800 


1.04 ± 0.33 


0.241 


1.00 ± 0.20 


0548-322 


Slew 


0.069 


0.161 






0.124 




6.289 


1.337 


9.558 


0.95 ± 0.05 






0605-085 


wp 


0.870 


1.547 





751 








0.546 


0.325 








0637-752 


wp 


0.651 


3.048 


1 


456 


0.200 


0.130 


8.500 


3.096 


2.854 


0.45 ± 0.48 






0716+714 


Uy 




0.933 





936 


0.165 




7.670 


2.395 


1.346 


1.77 ± 0.09 


0.139 


1.04 ± 0.33 


0735+178 


Uy,Slew 


>0.424 


1.590 


1 


085 


0.261 


0.143 


14.774 


1.783 


0.248 


1.34 ± 0.51 


0.368 


2.52 ± 0.83 


0736+017 


wp 


0.191 


1.593 


1 


041 


0.148 


0.076 


5.054 


1.462 


1.720 


1.82 ± 0.60 






0737+746 


Slew 


0.315 


0.019 












0.598 


1.083 


0.91 






0814+425 


Uy 


0.258 


2.615 





402 






1.531 


0.162 


0.041 


0.16 ± 0.78 






0820+225 


Uy 


0.951 


0.846 





133 






0.381 


0.070 


0.052 


1.05 ± 0.47 






0823+033 


Uy 


0.506 


0.976 





916 






3.958 


0.869 











•ntinued 



0828+493 


Uy 


0.548 


0.665 


0.409 






1.271 


0834-201 


wp 


2.752 


0.669 


0.211 








0836+710 


wp 


2.170 


0.578 


0.241 








0851+202 


Uy,Slew 


0.306 


2.173 


2.394 


0.791 


0.336 


22.177 


0859-140 


wp 


1.327 


0.767 


0.367 






1.377 


0906+430 


wp 


0.670 


0.924 


0.337 






1.306 


0923+392 


wp 


0.699 


4.067 


1.123 






1.905 


0927+500 


Slew 


0.188 


0.016 










0950+445 


Slew 


0.207 ? 


0.003 








0.016 


0954+556 


wp 


0.909 


0.936 


0.329 








0954+658 


Uy 


0.367 


1.068 


0.469 




0.109 


7.598 


1011+496 


Slew 


0.200 


0.247 








3.927 


1028+511 


Slew 


0.239 


0.037 










1055+018 


wp 


0.888 


1.519 


1.301 






5.055 


1101+384 


Slew 


0.031 


0.705 


0.323 


0.178 


0.134 


50.558 


1101-232 


Slew 


0.186 


0.058 








2.804 


1118+424 


Slew 


0.124 


0.032 










1127-145 


wp 


1.187 


2.376 


0.506 








1133+704 


Slew 


0.046 


0.264 




0.285 




9.241 


1144-379 


Uy 


1.048 


1.118 


0.946 


0.295 


0.116 


8.986 


1147+245 


Uy 




0.833 


0.380 


0.168 




10.205 


1148-001 


wp 


1.982 


0.459 


0.011 








1212+078 


Slew 


0.136 


0.085 










1215+303 


Slew 


0.237 


0.375 






0.158 


12.181 


1218+304 


Slew 


0.130 


0.051 








4.352 


1219+285 


Slew 


0.102 


0.908 


0.159 






11.124 


1226+023 


wp 


0.158 


35.419 


10.665 


2.204 


0.893 


88.981 


1239+069 


Slew 


0.150 


0.009 










1248-296 


Slew 


0.370 


0.004 










1253-055 


wp 


0.536 


8.557 


7.080 


0.240 


0.290 


6.041 


1255+244 


Slew 


0.141 


0.007 










1308+326 


Uy,Slew 


0.997 


0.761 


0.525 


0.425 


0.206 


7.857 


1320+084 


Slew 




0.010 










1332-295 


Slew 


0.513 


0.008 










1402+042 


Slew 


0.344 ? 


0.016 


0.018 






2.217 


1418+546 


Uy 


0.152 


0.946 


0.917 


0.217 


0.108 


12.737 


1421+582 


Slew 




0.006 










1424-418 


wp 


1.522 


0.940 


0.495 






8.651 



0.138 


0.035 





68 


± 





63 






0.152 


0.113 
















1.938 


0.819 





42 


± 





04 


0.315 


1.41 ± 0.18 


2.639 


1.063 


1 


50 


± 





17 


0.187 




1.405 


0.487 
















0.210 


0.088 





57 


± 





10 


0.275 




0.569 


0.712 


1 


26 


± 





13 






0.433 


1.905 
















0.113 


0.294 


1 


76 


± 





25 






0.405 


0.112 


1 


17 


± 





14 


0.020 


0.63 ± 0.23 


1.835 


0.158 





96 


± 


1 


31 


0.064 


0.85 ± 0.24 


1.526 


2.351 


1 


49 


± 





08 






0.743 


2.802 


1 


44 


± 





05 






0.558 


0.201 





93 


± 





44 






16.796 


37.333 


2 


10 


± 





08 


0.054 


0.58 ± 0.22 


0.601 


10.252 
















0.542 


3.896 
















1.060 


0.659 












0.794 


1.15 ± 0.36 


4.084 


2.681 


1 


51 


± 





10 






0.926 


0.663 


1 


67 


± 





38 






1.307 


0.047 


1 


18 


± 





79 






0.732 


















1.269 


0.747 
















2.899 


0.591 


1 


88 


± 





08 






1.337 


10.324 


1 


22 


± 





03 






3.047 


0.409 


1 


24 


± 





16 


0.010 


0.27 ± 0.39 


29.834 


12.074 





89 


± 





05 


0.234 


1.40 ± 0.10 


0.057 


0.806 
















0.096 


2.376 
















2.073 


1.246 





83 


± 





04 


1.612 


0.89 ± 0.06 


2.420 


10.184 
















1.818 


0.134 


1 


04 


± 





33 






0.114 


5.016 
















0.076 


0.360 


1 


14 


± 





21 






0.571 


0.874 


1 


85 


± 





17 






1.970 


0.306 


1 


14 


± 





19 






0.142 


1.342 
















1.050 














5.50 


1.6 ± 0.4 



ible 1 - continued 



1440+122 


Slew 





162 


0.037 












1504-166 


wp 





876 


0.874 





267 








1508-055 


wp 


1 


185 


0.880 





153 






1.333 


1510-089 


wp 





361 


2.063 


1 


942 






5.173 


1514-241 


Uy 





049 


1.907 


1 


084 


0.233 


0.167 


23.579 


1517+656 


Slew 






0.033 












1519-273 


Uy 






1.958 





336 






1.049 


1533+535 


Slew 






0.008 












1538+149 


Uy 





605 


1.452 





319 


0.080 




2.485 


1544+820 


Slew 






0.036 












1553+113 


Slew 





360 


0.497 






0.208 


0.099 


24.440 


1610-771 


wp 


1 


710 


1.073 





298 








1611+343 


wp 


1 


404 


0.854 





469 






0.978 


1633+382 


wp 


1 


814 


0.763 





461 






1.606 


1641+399 


wp 





594 


4.266 


2 


927 


0.724 


0.281 


7.825 


1652+398 


Uy,Slew 





034 


1.383 





313 


0.105 


0.066 


41.250 


1727+502 


Slew 





055 


0.152 










3.737 


1741+196 


Slew 





083 


0.209 












1741-038 


wp 


1 


046 


1.455 





806 








1749+096 


Uy 





322 


1.384 


1 


646 


0.214 


0.078 


5.031 


1749+701 


Uy 





770 


1.023 





217 








1803+784 


Uy 





684 


1.556 





674 


0.307 


0.143 


7.511 


1807+698 


Uy,Slew 





051 


1.646 





822 


0.273 


0.154 


27.842 


1823+568 


Uy 





664 


0.998 





559 






4.117 


1853+671 


Slew 





212 


0.010 












1928+738 


wp 





360 


2.239 





622 








1954-388 


wp 





630 


1.091 


1 


179 






5.950 


1959+650 


Slew 





048 


0.242 












2005-489 


Uy 





071 


1.148 





515 


0.260 


0.237 


30.090 


2007+777 


Uy 





342 


0.939 





630 








2052-474 


wp 


1 


489 


0.770 





293 








2106-413 


wp 


1 


055 


0.921 





578 








2128-123 


wp 





501 


1.233 





413 


0.080 




3.280 


2131-021 


Uy 





557 


1.362 





385 






1.029 


2134+004 


wp 


1 


936 


3.020 





331 






1.353 


2145+067 


wp 





990 


1.844 


1 


734 








2155-304 


Slew 





117 


0.284 





241 


0.088 


0.096 


37.697 


2200+420 


Uy 





069 


4.462 


1 


959 


0.467 


0.239 


37.519 



0.482 


1.617 












0.337 


0.203 












0.795 














1.276 


0.718 


0.90 


± 


0.16 


0.269 


1.51 ± 0.36 


5.474 


0.607 












1.445 


3.410 












0.217 


0.417 


1.37 


± 


1.78 






0.307 


1.457 












0.485 


0.092 


1.05 


± 


0.90 






0.640 


2.596 












6.978 


7.656 












0.535 














0.553 


0.194 


0.76 


± 


0.06 


0.317 


1.00 ± 0.20 


0.401 


0.258 


0.53 


± 


0.08 


0.739 


1.03 ± 0.09 


1.135 


0.914 


0.85 


± 


0.23 






11.040 


8.475 


1.63 


± 


0.05 






1.668 


3.707 


1.39 


± 


0.08 






0.845 


3.079 












0.753 


2.213 








0.921 


2.00 ± 0.40 


1.200 


0.129 


0.45 


± 


1.43 






1.079 


0.193 


1.44 


± 


0.71 






0.683 


0.324 


1.42 


± 


0.45 






5.099 


0.317 


0.81 


± 


0.30 






0.198 


0.241 


0.44 


± 


0.38 






1.046 


0.801 












1.308 


1.295 


1.33 


± 


0.19 






1.113 


0.353 












17.220 


9.801 












12.651 


5.333 


1.94 


± 


0.06 


0.117 




0.819 


0.158 


0.75 


± 


0.56 






0.135 


0.433 








0.296 


1.40 ± 0.40 


0.053 


0.312 












3.167 


1.917 












0.145 


0.051 


1.05 


± 


0.46 






1.268 


0.214 


0.82 


± 


0.42 






1.417 


0.392 












18.269 


45.583 


1.34 


± 


0.03 


0.124 


0.71 ± 0.24 


5.527 


0.936 


1.92 


± 


2.46 


0.338 


1.20 ± 0.30 



ible 1 - continued 



2203-188 


wp 


0.618 


2.290 


0.210 






0.546 


0.184 












2223-052 


wp 


1.404 


1.445 


3.346 


0.841 


0.347 


10.821 


1.615 


0.292 


1.09 ± 0.23 








2230+114 


wp 


1.037 


1.463 


0.898 


0.192 




2.338 


0.701 


0.486 




0.437 


1.60 ± 


0.20 


2240-260 


Uy 


0.774 


0.581 


0.129 






1.890 


0.295 


0.062 


0.79 ± 0.40 








2243-123 


wp 


0.630 


1.298 


0.776 






0.837 


1.997 


0.963 










2245-328 


wp 


2.268 


0.397 


0.097 








0.329 












2251+158 


wp 


0.859 


3.912 


4.057 


0.188 


0.115 


10.285 


1.407 


1.082 


0.62 ± 0.04 


0.699 


1.18 ± 


0.08 


2254+074 


Uy 


0.190 


0.975 


0.249 


0.155 




4.966 


0.619 


0.105 


1.89 ± 0.61 








2321+419 


Slew 


0.059 


0.018 










0.441 


0.732 










2326-477 


wp 


1.302 


0.856 


0.415 








1.290 


0.245 










2343-151 


Slew 


0.226 


0.007 










0.073 


0.853 










2345-167 


wp 


0.576 


2.032 


0.916 


0.093 




1.156 


0.314 


0.216 











eferences of Table 1: 

1 the references on the data collected for the 126 sources are reported in the Notes to Table 5. 
otes to Table 1: 

' source detected only by Whipple and HEGRA, in the TeV band. 



